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£h| Abstract 
CD 

Motivated by the considerable success of alternative theories of gravity, we consider the toy 
model of a higher derivative Lagrangian theory, namely the Pais-Uhlenbeck oscillator studied 
C\) ■ in a recent paper by Hawking & Hertog. Its Euclidean Path Integral is studied with a certain 

detail and a pedagogical derivation of the propagator, which makes use of a Theorem due to 
Forman, is consequently proposed. 

1 Introduction 

t> 
O 
ON 

o 



The discovery made by Rless and Perlmutter and respective collaborators [1,2] that the universe 
where we live is expanding with an accelerating rate, probably represents the greatest challenge 
of the century that Nature has provided to theoretical physicists (see, for example [3] ) . Among 
the large number of models and physical mechanisms proposed in order to explain the accel- 
erating era of the universe, are of some interest the so-called /(i?)-theories of gravity (see, for 
example [4]). Loosely speaking, they represent alternative theories of gravity where, in place of 
the Einstein-Hilbcrt action 

Ieh = ieb J dx ^R 

the Ricci scalar is substituted by some appropriately chosen function /(i?) (or other higher 
order invariants like the Gauss-Bonnet invariant as in [5-9]). ^From this perspective, the present 
acceleration era is a manifestation of the new, more involved, geometry of the universe. This 
represents a conceptual difference from models where the acceleration is driven instead by the 
presence of exotic fields (quintessence, phantom field) which dominate the matter content of 
spacetimc. 
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We should point out that the idea of introducing a correction to the Einstein-Hilbert action 
in the form of f(R) = R + R 2 was proposed long time ago by Starobinsky [10] in order to solve 
many of the problems left open by the so-called Hot Universe Scenario 1 . 

Nowadays, /(i?)-theories of gravity are understood as toy models without any intention of 
dcfinitcncss, that is, useful playgrounds where new physics, possibly related to observable events 
of our universe, can appear. Their interest grew up in cosmology with the appearance of the 
papers [11, 12]. 

From a mathematical point of view, /(i?)— theories of gravity, polynomial in the scalar 
curvature and with degree higher than one of course (as the Starobinsky model), are known 
to be asymptotically free [13] and renormalizable [14-16]. However, the introduction of higher 
derivatives of the metric seems to lead to ghosts, states with negative norm which we think are 
able to spoil any QFT from physical interest. A recent paper by Hawking & Hertog [17] has 
revitalized the interest in such kind of theories since, as shown by the authors, starting from a 
positive definite action, one can give meaning to the Euclidean path integral as a set of rules 
for calculating probabilities for observations. With reference to the specific model of a fourth 
order Lagrangian, namely that of the Pais-Uhlenbeck oscillator (PU) [18] , it is shown that, 
paying the prize of losing unitarity, one can never produce or observe negative norm states. 
This result proves the goodness of the model and justify in our opinion the attention received 
recently in [19]. It is clear that the first step toward the comprehension of the PU quantum 
mechanical oscillator is represented by proper evaluation of its propagator. Both [17, 19] have 
performed such calculation. Our results, modulo a normalization constant, are in agreement 
with the ones presented in [19], also if different techniques have been implemented. Aim of the 
present paper is to evaluate explicitly step by step the propagator of the PU oscillator. 



2 Path Integral Representation of PU Oscillator Propa- 
gator 

Let us consider the one-dimensional PU Lagrangian [18], tel, 

L PU = U d 4) -V{ q )-^- , with m,a>0. (2.1) 



2\dt J V1/ 2 V* 2 J 
Its Euclidean version, obtained by Wick rotating the time t, i.e. t — > —it, is 

- £ -K£)' + ™ + t (*£)'• 

For brevity, we shall denote with an overdot the derivative with respect to r. 
Set h = 1 , we can formally write the propagator as 

Bqe l ^dtL PU w^k z T ( A ) := f Bqe -i E [ q ] (2.3) 

A J A 

with the Euclidean action given by 

IeM = £ dr Qg 2 + V(q) + ^-f^j . (2.4) 

The Euclidean action turns out to be positive definite as far as V(q) is non negative and the 
propagator (2.3) as explained in [17] can be used, at least in the Gaussian approximation, 



1 This in turn had the consequence of introducing an accelerating expansion in the primordial universe, so that 
the Starobinsky model can be considered as the first inflationary model. 
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to extract probabilities for physical observations. Here, Dq represents the formal functional 
measure and A, the boundary conditions necessary to give a meaning to a formal path integral. 
As well known, in the usual second order theory, it would be sufficient to specify q on the initial 
and final time slice in order to make the propagator well defined. However, the present theory 
is of order higher than two so that extra boundary conditions are expected to be involved in the 
definition of (2.3). As already proposed in [20], and fully explained in [17], the right choice is 
provided by 



A: 



q(0) = q , q{T) = q T 



q(0) = q , q(T) = q T , (2.5) 



since, any condition on q would make otherwise the action infinite. 

Established the boundary conditions, the propagator (2.3) can be re-written in the more 
evocative form 



Z T (A) = (q T ,qT-,T = T\q ,q Q :T = 0). 



(2.6) 



Now, one may formally procedes by splitting q into a "classical" part, q c i, and a quantum 
fluctuation q, i.e. 



(2.7) 



q c i solves the classical EOM obtained by SIe = with boundary conditions (2.5). From (2.7) and 
(2.5), it turns out that quantum fluctuations have to satisfy the following boundary conditions, 



A: 

The Euclidean action (2.4) becomes: 



3(0) = = q(T) 



q(0) = = q(T) 



(2.8) 



lE[q c l +q] 



lE[q c i] + Ie 



[q c i,q\ + 

Jo 



dr (qq c i + m 2 qq c[ + a 2 qq c[ 



pi 



lE[q c i] +I E [qci,q] + drq 
Jo 



( 





q c iq- 



! d 4 q cl _ d 2 q ci 



a 2 q c iq 



o? q ciq 



dr 2 

T 




+ V'(q cl ) 



+ 



where 



AlclA] = 

Jo 



dr 



\k 2 + ^'q 2 + W{q cU q) 



(2.9) 



(2.10) 



with 



W(q cl ,q) = V(q cl + q) - qV'(q cl ) - V(q cl ) . 



(2.11) 



Notice that, by extremizing the action, we get the EOM for the classical solution, namely. 



EOM: 



a 



2 d 4 q cl d 2 q cl 



+ V'(q cl )=0. 



(2.12) 



dr 4 dr 2 

In (2.9) the integral vanishes on-shell, while the boundary term vanishes upon using (2.8). To 
go on, we make the usual Gaussian approximation or cquivalently, we may restict to quadratic 
potentials of the type V(q) = ,J irq 2 - Thus, W(q c i,q) = W(q) = ,J irq 2 , with m 2 constant. The 
Euclidean action after the splitting (2.7) and in reason of this consideration neatly separates 
into 



IeM 



iE[q c i] + Ie[q] 



( = 8) lE ^ + \j n dT $ 



(r) 



d 4 



d 2 



a 



dr 4 dr 2 



+ 77T 



q(r) 



(2.13) 
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and the classical action has been explicitly evaluated in [17]. 

The propagator assumes the form of a Gaussian integral in the quantum fluctuation variables: 



(2 13) 

(Qt,C[t;t = T\q ,qo]T = 0) = cxp-I E {q cl ] 



od A d 2 , 

a d^~d^ +m 



q{r) (2.14) 



As a consequence, one has to give a meaning to the formal Gaussian path integral. To this 
aim, let us denote 



2 d A d 2 2 
Kx:=a d^-d^ +m 



(2.15) 



as the fourth-order differential operator in L2(0,T), defined in the dense domain D(K) := 
{f,K x f e L 2 {0,T) | X(f) = 0} with suitable boundary conditions X. 

Let us try to determine X such that K is a self-adjoint operator in Z/2(0,T). To this purpose, 
suppose /(") (n = 0, . . . , 4) be in D{K): 

d A d 2 



[ dTg{r) 
Jo 



2 u. u, 2 



/ 



ff(r)/(r) + 



dr 4 rfr 2 

+.9/ + a 2 (g /' + ^/ -■§/) 
= (Kg,f) + B(0,T). 
K is certainly symmetric if and only if 

B(0,T) =0. 



(2.16) 



(2.17) 



This condition is realized, among the others, by the functions /^( T )i 3^( T ) ( n = 0>---)4) 
absolute continuous in [0,T] s.t. 

/(0) = = /(T) & /(0) = = f(T) ■ 3 (0) = = g(T) & . 9 (0) = = g{T) (2.18) 



or 



/(0) = = /(T) & /(0) = = f(T) ; 3 (0) = = g(T) & 3(0) = = g(T) . (2.19) 

Thus, we may invoke the general Von-Neuman-Krein method to find all the self-adjoint exten- 
sions of a symmetric operator. In our case, one has to find the L 2 (0, T) solutions to equation 



K^u±(t) — ±iu±(r) . 



(2.20) 



Since (0,T) is compact, it turns out that the defect indices are (n + ,n_) = (4,4) meaning 
that all self-adjoint extensions are parametrized by 4 x 4 unitary matrices. For our pourposes, 
we simply observe that physical boundary conditions (2.18), corresponding to (2.5) - (2.8), 
can be represented by an unitary matrix, and the operator K is self-adjoint. This is also 
confirmed by the fact that the associated spectral problem is well defined, eventhough the 
equation which implicitly defines the eingenvalues is highly trasccndcntal. In fact, for example, 
even the spectrum of the simplest fouth order operator := a 2 -^ with boundary conditions 
(2.18) , is analytically unaccessible due to the impossibility of solving the trascendent equation 2 



2 Said A„ the n-th eigenvalue of K^°\ then z„(A n ) := i/aA^' 



1/4 



4 



cos(Tz„) = scch(Tz„). 

On the other hand, it is also easy to show that the boundary conditions (2.19) defines also a 
self-adjoint extension. But in this case, the spectral problem is much more easier to handle. In 
fact, let us denote by K 

K = L a + m 2 + a 2 Ll , with (2.19) as BC , (2.21) 
where L = — -§pz- The eigenfunctions are sin (^) and the the spectrum is 

a(K) := {A„ := (^) 2 + m 2 + c? {^f) , n e N} . (2.22) 



The problem is that boundary conditions (2.19) are "unphysical" , in the sense that they do 
not enter in the expression of our propagator. Nevertherless, as we will see, they will play an 
important role. 

We conclude this Section with the final form of PU propagator, obtained performing the 
Gaussian integral (2.14): 



/ 2ir 

(q T , qr\T = T\ q , q ; r = 0) = -y/ pc cxp {-I E [q c l]) ■ (2.23) 

3 Regularization of Functional Determinants 

The goal is to give a rigourous meaning to the formal functional determinant which appears 
in (2.23). This is also called prefactor. In [17], the authors have performed the calculation of 
the prefactor under the (implicit) assumption that the Van-Vleck-Pauli method is valid even 
for the higher-order dynamical system. This is not obvious, since the Pauli theorem is based 
on Hamilton- Jacobi theory for ondinary second order Lagrangian. In our case, the Hamiltonian 
formalism is quite different from the ordinary one, known as Ostrogadsky formulation. A second 
attempt has been recently proposed by [19] but we reserve to comment it later. 



In our approach, we shall regularize the functional determinant by zeta-function regularization 
[21-25]. A regularization is necessary since functional determinants arc formally divergent. We 
recall that a simple regularization for the functional determinant associated with an elliptic 
operator L may be chosen as 



In Det L(e) 



-f 

Jo 



dt 



1 



r(i + e) 

where the zeta-function is defined by means of the Mellin-like transform 
1 



C(s\L) 



T(s) 



dt t s - lr Tre- tL 



& 



C 



(3.1) 



(3.2) 



For a Q order differential operator in D-dimensions, the integral is convergent in the region 
Res > |j where the function C( s l^) is analytic. It is possible to show that C( s l-^) can be 
analytically continued in the whole complex plane and it is regular at s = 0. Thus, by expanding 
(3.1) in Taylor's series, we obtain 



In DetL(e) 



4c 



(3.3) 
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and the regularised functional determinant associated with L can be defined by taking the finite 
part in the limit e — > 0, that is 



In Det L = -(' (° 1 7^) = -C'(0|i) - In V? C(0|i) , 



(3.4) 



leading to the usual zeta- function regularisation prescription [21,22]. The divergences are gov- 
erned by the computable coefficient C(0|L), which does not depend on the arbitrary scale pa- 
rameter [i. 



Within this direct zeta- function regularization, in order to calculate Det-KT^, we need to know 
explicitly the spectrum. As alreday stressed, often one does not know explicitly the spectrum. 
Thus, we are forced to make use of a powerful theorem proved by Forman in [27], which is 
a generalization of Gelfand-Yaglom [28] and Levit-Smilansky theorems Levit-Smilansky theo- 
rems [29]. Adapted to the case at end, we may state it in the following way 3 : 



Theorem (Forman, 1987): Let K, K be operators defined in [0,T] of the form 
d n / d n ~ 2 



dr n ' \dr n - 2 



K = P ^ + {d^)- ^ 
Let A be any boundary condition represented by matrices (M,N), s.t. 

MO) 

Ml ••• |+iV| 1=0. (3.6) 



(3.7) 




^""^(O) 
Then, 

BetK A _ det(M + NY K {T)) 
V>ctK A ~ det(M + NY R {T)) 



where (M,N) are matrices defining other boundary conditions smoothly related to (M,N). 
Above, for any h such that Kh — 0, Yk{t) acts as 





Y K (r) ■■■ . (3.8) 



A solution of this equation is given by 



Yk(t) 



( u i( T ) 
ui(r) 
iii(r) 



u 2 (t) 
m 2 (t) 
1*2 (r) 
u 2 {t) 



«3(t) 



U 4 (t) \ 

u 4 (t) 

1*4 (t) 

u 4 (t) J 



(3.9) 



3 What we report here is a simplified version of Proposition 3.9 in [27]. 
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where Kuj = for j = 0, . . . 4, satisfying the initial conditions 



«i(0) = 


1 , 


Uj (0) 


= 0, 




n 2 (0) = 


1 , 


uj(0) 


= 0, 


2 , 


« 3 (0) = 


1 , 


Uj(0) 


= 0, 


.7^3, 


■« 4 (0) = 


1 , 


«i(0) 


= 0, 


i ^ 4. 



The role of K which appears in the Theorem is played by the operator (2.15). The boundary 
conditions (2.8) can be put in correspondence with matrices 4 

«=f'!;> SI. « = f.° SV (3") 



o o ; ' "V % , 

Finally, in our case, the matrix Yk(T) can be easily computed. The general solution of Ku = 
reads 

Uj(r) = Aj sinh(Air) + Bj cosh(Air) + Cj sinh(A 2 T) + Dj sinh(A 2 T) , (3.12) 

with 



Al ' 2 - V 2^ ' (3J3) 

In the oscillatory regime, that is 2am < 1, the roots Ai j2 are real as can be checked expanding 

the double radicals above. 

Imposing the initial conditions (3.10), one has 

A 2 A 2 
u i( r ) = x2 2 12 cosh(Air) - 1 cosh(A 2 r) , 



A 2 A 2 

u 2(t) = Ai( - A | 2 _ A 2^ sinner) - ^ sinh(A 2 r) 
^"^2 cosh(Air) + 

A 2 A^ A 2 — A 1 

A^A^^^ + W-A 2 ) 
Therefore, the right hand side numerator of (3.7) is 



Mt) = -T2 72 COsh(A!T) + -3 -j CO S h(A 2 T) 

\ 2 — A 1 A 2 — A : 

u 4 (r) = - ^^2 ^ 2 ^ sinh(Air) + / 2 sinh(A 2 T). (3.14) 



det(M + a^(t)) (3 = 1} U3 (r)u 4 (r) - U4 (r)w 3 (r) 



2-2 cosh(AiT) cosh(A 2 T) + 1 2 sinh(AiT) sinh(A 2 T) 



(A 2 - A 2 ) 2 V ' ' K ' AiA 



2 



i.-i.i.'i) a 3 [ 1_ 2 / yi + 2amT \ _ 1 2 / -y/1 - 2amT \ 

1 + 2am bm V 2a / 1 - 2am &m \ ^ / 



m 



.(3.15) 



In order to compute completely the propagator, we need to find the fourth-order differential 
operator which will play the role of K in Forman's theorem. The requirement on the candidate 
is that we must be able to compute its funtional determinant in an independent way. Loosely 



4 The choice is highly non unique of course. The final result, however, seems not to depend on this large freedom. 
Cfr. [27]. 
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speaking, there are two ways to compute the functional determinant of an operator, that is: (i) 
taking the product of its eigenvalues (if you know the spectrum) ; (ii) using some "smart" math- 
ematical theorem. We have already decided to follow the latter to compute Det K, something 
which forces us to take the former way for Det K. 

As we have seen, Forman's theorem lets us to play with the boundary conditions both of K and 

K. This suggests us to choose K formally equal to if, with unphysical boundary conditions 

(2.19). We recall that in this case we known the spectrum, given in (2.22). 

A quick way to compute the functional determinant Det K may be the following. Zeta-funtion 

regularized determinants suffer, in general of the so-called multiplicative anomaly [30-32], 

namely 

Det( AB) = Det A Det B e a< - A ' B ^ . (3.16) 
The quantitya(A, B) is called multiplicative anomaly. In our case, by noticing that 



- / T 1 - VI - 4o< 2 m 2 \ / T 1 + Vl - 4a 2 m 2 \ - - 
K=laL Q + Ul + \=K-K+, (3.17) 

and observing that, since K is an ordinary fourth order differential operator, the multiplicative 
anomaly is vanishing (cf. [32]), we have 

Detif = (Det if_) (Det if +). (3. 18) 

In this way, K± are just second order differential operators obtained by shifting the simple Lq 
with a constant term. The calculation of (Det_K±) is standard and we are not going to repeat 
it here. The result is 



Det if 



mT 2 



2 I Vl + 2amT\ 2 f VI ~ 2amT\ 

-nili" ! i - sinn 



J- Smh I Ya ) 



(3.19) 



We may confirm this result (and in turn, the absence of multiplicative anomaly) by evaluating 
directly the determinant of K. This target can be accomplished by implementing a standard 
trick which consists in the observation that the derivative with respect to m 2 of ln(Dct K) is a 
well defined quantity, namely 

'' ln(Detif) = - Jz ^tr(lnif) 

OO 



d(m 2 ) v ; d(m 2 ) 

d 

d^ 2 ) 

OO 1 

^ K 



n=l 

< oo (3.20) 



n=l 

Thus, the problem has been reduced to compute the convergent series 



71=1 

1 



E 



„=i K -) 2 + m 2 + a 2 (?f) 4 



T2 f 1 (3 21 ) 



Vl - 4a 2 m 2 \ ^ (7rn) 2 + z 2 _ ^ {im) 2 + z\ 

v \n—l v ' n—1 x / * , 
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where for simplicity we have introduced the quantities 



^ /l±Vl-4a 2 m 2 T , /- , N 

z± :=T\ - , = — Ul + 2am ± Vl - 2am) . 

V 2a z 2a 



In terms of these new quantities, 

-{2z±)dz±. 



,, 2 , Vl — 4a 2 m 2 
a(m J = T- 



y2 

Recalling the Mittag-Lefflcr expansion of the coth function, i.e. 

1 °° 1 
cothz = - + 2z\^ - — — -, 

z ^ im 2 + z 2 

n=l K ' 

it turns out that 

ln(Det^) = J dz- ^coth(z_) - — ^ + J dz + ^coth(z+) - — ^ 

/sinh(z + ) sinh(z_) \ 
n V z + z - ) ' 
so that the final result coincides indeed with (3.19). 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



Now that we know Det K from independent considerations, we can close the circle just computing 
the denominator of (3.7). Since K and K are formally the same, the Y matrix does not change. 
What changes are the boundary conditions which now are represented, for example, by 



M = 



( 1 o 





\ 



\ 

-1 



j 



N = 



( 








\ 














1 











V o 





1 





(3.26) 



The right hand side denominator of (3.7) becomes 
det(M + NY R (T)) = u 2 (T)u 4 (T) - u A (T)u 2 (T) 
■ sinh(AiT) sinh(A 2 T) 



(3.14) 1 



A1A2 
(3.13) a 
m 



sinh 2 r^^Vsinh 2 
2a 



VI - 2amT\ 
2a ) 



l - , , ■ (3-27) 

\ Za J \ la J 

Inserting (3.27), (3.15) and (3.19) into (3.7), the functional determinant of the K operator with 
boundary conditions (2.5) - (2.8) is 



DetK A 



a 



mT 2 



1 + 2am 



sinh 2 f VT ™^ + 



2a 



1 



2am 



■ sinh' 



/ VI -; ■'<"•'<' r 



2a 



(3.28) 



This is the main result of the paper. A dimensional analysis tells us that, in natural units, 
[a] = [T] = (mass) -1 , so that [Det if] = (mass) -2 . Up to a minor adjustment of the numeric 
factor, (3.28) is the prefactor presented in [19], equation (9). 
It is straightforaward to check that, in the small T limit, 



rp2 rpi 

Det K a ks— + ? + 0(T 6 

A 12 180a 2 V 



T< 1 , 



(3.29) 
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but for t € [0, oo), 
Det K A = 



a 
m 



1 exp (VI + 2am^) 



1 + 2am 



(2T/a) 2 
1 



+ 



exp (\/l — 2am^) 



1 — 2am 



(2T/a)* 



T > 1 , (3.30) 



something which shows that the ground state probability amplitude, being proportional to 
Det-ftTT 1 , is indeed exponentially suppressed and makes the Euclidean theory well defined. 
Finally, notice that the vanishing "mass" term case is 

v 4 



DetK A 



>0 OL 



2 1- cosh 



T 

a 
(i) 



T ■ , 
H sinh 

a 



=: Det if 



(i) 



(3.31) 



where we have denned the operator K A ' := a 



dr 4 



d^ 
dp 



with boundary conditions (2.5) - (2.8). 



a 2 ^ in t e [0,T] (operator 



Given (3.31), one can evaluate, for example, the determinant of the most elementary fourth- 
order differential operator of physical interest, namely := 
whose spectrum, as said above, is unaccessible to us), which is 

Det K {0) - — 



A 12 ' 

We leave the details of the computation to the interested reader. 



(3.32) 



4 Conclusions 

To summarize, we have devoted the paper to the computation of the quantum mechanical Eu- 
clidean propagator of the one-dimensional PU oscillator. As we have seen, the PU Lagrangian 
gives rise to a fourth order differential operator whose functional determiant, Det K A , has repre- 
sented the main difficulty along the way. As showed by Hawking & Hertog, the correct boundary 
conditions A we have to impose in order to give meaning to the Euclidean propagator, are pro- 
vided by equation (2.5). By making use of Forman's theorem, we have been able to evaluate 
DetK A , so that the main result of the paper can be stated in the following way, 



2ir 

(q T , q T ;r = T\ q , q ;r = 0) = \j p^r- exp {-I E [q c l\) , (4.1) 



where 

Det K A = 



a 3 1 . , o /VI + 2amT\ 

■ sinh + 



mT 2 



1 + 2am V 2a 
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1 . , 2 ( Vl — 2amT\ 
■ sinh 



1 — 2am \ 2a J 

and the classical Euclidean action can be found in [17], equations (A5)-(A6). 



(4.2) 



As a byproduct, we have shown how Forman's theorem can be usefully implemented to give other 
functional determinants of potential interest. To this regard, Forman's theorem proves to be 
an instrument of extraordinary power. The generalization to higher order quadratic Lagrangian 
L(q,q^\ - ■ ■ ,q^) is immediate, eventhough the computations become of course much more 
involved. 

We conclude with the observation that the method described in this paper may also be useful 
in the so-called Horava-Lifshitz non-relativistic renormalizzable theory of gravity [33], where 
higher spatial derivatives appear indeed in the Lagrangian of the theory and in inflationary 
cosmology [34]. 
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